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ABSTRACT 


Following the Challenger tragedy, an evaluation of the integrated main propulsion system 
flight data revealed a premature decay in the hydrogen external tank ullage pressure. A 
reconstruction of predicted ullage pressure versus time, to explore this anomaly, indicated an 
inconsistency between predicted and measured ullage pressure starting at approximately 65.5 
seconds into the flight and reaching a maximum value between 72 and 72.9 seconds. This 
discrepancy could have been caused by a hydrogen gas leak or by a liquid hydrogen leak that 
occurred either in the pressurization system or in the external tank. The corresponding leak rates 
over the time interval from 65.5 to 72.9 seconds were estimated to range from 0.28 kg/s (0.62 lbm/s) 
± 41 percent to between 0.43 and 0.51 kg/s (0.94 and 1.12 lbm/s) ± 1 percent for a gas leak and from 
72.9 kg/s (160.5 lbm/s) ± 41 percent to between 111.6 and 133.2 kg/s (245.8 and 293.3 lbm/s) ± 1 
percent for a liquid leak. No speculation is made to ascertain whether the leak is liquid or gas, as 
this cannot be determined from the analysis performed. Four structural failures in the hydrogen 
external tank were considered to explain the leak rates. A break in the 5-centimeter (2 inch) 
pressurization line, in the 13-centimeter (5 inch) vent line, or in the 43-centimeter (17 inch) feedline 
is not likely. A break in the 10-centimeter (4 inch) recirculation line with a larger structural failure 
occurring in the 72- to 73-second time period, the time of the visibly identified premature pressure 
decay, does seem plausible and the most likely of the four modes considered. These modes are not 
all-inclusive and do not preclude the possibility of a leak elsewhere in the tank. 


INTRODUCTION 


Subsequent to the Space Transportation System (STS) 51-L accident on January 28, 1986, the 
Propulsion and Power Division at the NASA Lyndon B. Johnson Space Center (JSC) formed an 
investigative team to assess and evaluate the main propulsion system (MPS) performance. This 
team concentrated on the integrated MPS consisting of 

1 . The Orbiter main propulsion system 

2. The Space Shuttle main engines 

3. The external tank 

Flight data were scrutinized for anomalies that occurred during the Challenger flight. 

During the investigation, a review of the MPS flight data showed that a premature decay in 
the hydrogen external tank ullage pressure started seconds before the catastrophe. Under normal 
conditions, once the hydrogen external tank pressurization system stabilizes after the initial 40 
seconds of flight, the system operates such that the ullage pressure in the external tank is 
maintained within the bandwidth shown in figure 1. 
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Normal 

bandwidth 


Figure 1.- Hydrogen external tank ullage pressure, STS 51-L. 

The pressure decay indicated in figure 1 could be indicative of a hydrogen leak in the external 
tank. The accuracy of this postulation needed to be determined because the decay could be 
explained also if a pressurization system malfunction had occurred. For example, if a valve that 
permits pressurant gas to flow to the tank was never commanded open, an insufficient quantity of 
pressurization gas would be recirculated to the tank and a subsequent drop in ullage pressure would 
result. 


In compliance with the NASA’s publication policy, the original units of measure have been 
converted to the equivalent value in the Systeme International d’Unites (SI). As an aid to the 
reader, the SI units are written first and the original units are written parenthetically thereafter. 

The author wishes to extend appreciation to Regina Rieves for her technical support and 
advice, and to Dr. Gene Ungar and Dr. A1 Feiveson for the technical inputs during the review of this 
report. The author also wishes to extend thanks to Warren Brasher for his supervisory support of 
this project. 


SUMMARY OF RESULTS 


The STS 51-L hydrogen external tank ullage pressure was reconstructed using flight data 
(app. A) independent of measured pressure such that a comparison of predicted to measured ullage 
pressure would either verify or discount the presence of a leak. If no leak was present, then the 
predicted ullage pressure, which uses actual valve position flight data (i.e., it was assumed that 
flight data accurately indicated the position of each valve), should match the measured ullage 
pressure. The pressures would match regardless of whether or not an additional valve that would be 
required to open to prevent the decay was actually commanded open. Thus, if predicted and 
measured ullage pressures do match, one could then attribute the premature decay to a malfunction 
of the main propulsion system and delve further into the control system dynamics of the 
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pressurization system. If a leak was present, however, the predicted ullage pressure would deviate 
from that measured. 

Ullage pressure was reconstructed using data from flights STS 61-A and STS 51-F as a 
baseline. A baseline was desired for two reasons. First, since the thermodynamics inside the ullage 
are complicated and undefined, data from flights STS 61-A and STS 51-F were needed to correlate 
the bulk ullage temperature, an important parameter required in the reconstruction. Second, 
consistent results between several flights would affirm the integrity of the analysis and confirm the 
validity of the results. Flights STS 61-A and STS 51-F were chosen because their throttling profiles 
(liquid hydrogen consumption rates) are similar to those of STS 51-L, and it was assumed that 
similar thermodynamics would be present. 

Results of ullage pressure reconstruction for flights STS 51-L, STS 61-A, and STS 51-F, from 
40 to 73 seconds mission elapsed time (MET), are shown in figures 2, 3, and 4, respectively. The 
timeframe prior to 40 seconds was not considered since the period of interest included only the final 
seconds of flight. As shown in figures 3 and 4, there is agreement to within ±1.7 kPa ( ± 0.25 psia) 
between predicted (reconstructed) and actual (measured) ullage pressure, from the start of recon- 
struction at 40 seconds to the end of reconstruction at 73 seconds, for both flights STS 61-A and STS 
51-F. These results are consistent with a no-leak assumption. On the other hand, as shown in 
figure 2, flight STS 51-L shows agreement to within ± 1.7 kPa (±0.25 psia) until approximately 64 
seconds. At this time, the predicted ullage pressure began to deviate significantly from actual 
ullage pressure, with the predicted ullage pressure of increasingly greater magnitude than the 
actual ullage pressure. The result verifies the postulation that a leak was the source of the 
premature ullage pressure decay. 

The deviation of the ullage pressure from the prediction can be explained by one of two 
phenomena. One, the ullage volume is increasing faster than measured, implying a liquid leak in 
the tank. Two, the amount of pressurization gas actually entering and remaining in the tank is less 
than measured, implying a gaseous hydrogen leak either in the pressurization line leading into the 
tank or in the tank itself. 

With an apparent leak identified and the leak initiation time estimated at 65.5 seconds, the 
next objective was to quantify the leak rates. Actual and predicted ullage pressure rise and decay 
rates, which are dependent upon hydrogen mass quantities in the tank, were chosen as the criteria 
from which leak rates could be quantified. If comparison of actual and predicted pressurization 
rates showed a deviation, this mismatch could be corrected by incorporating a leak into calculations 
of predicted pressure rise and decay rates. To calculate a hypothetical liquid leak, no gas leak is 
assumed. Likewise, to calculate a hypothetical gas leak, no liquid leak is assumed. 

The time period in question, from 62 to 73 seconds, was broken into three separate intervals, 
and actual and predicted pressurization rates were calculated for each of these intervals and 
compared. Leak rates were computed at average times on each interval and are shown in figures 5 
and 6. If a gas leak is assumed, calculations show an initial leak of 0.28 kg/s (0.62 lbm/s) ± 41 
percent at 65.5 seconds which increases to between 0.43 and 0.51 kg/s (0.94 and 1.12 lbm/s) ± 1 
percent in the 72- to 72.9-second time interval. If a liquid leak is assumed, calculations show an 
initial leak of 72.9 kg/s (160.5 lbm/s) ± 41 percent at 65.5 seconds which increases to between 111.6 
and 133.2 kg/s (245.8 and 293.3 lbm/s) ± 1 percent in the 72- to 72.9-second time interval. In both 
cases, the leak remains relatively constant from 65.5 to 72.0 seconds and increases significantly at 
72 to 73 seconds. This sharp increase in leak rates parallels the premature decay that was initially 
identified in the flight ullage pressure data of figure 1 . 

Finally, a matrix of several possible structural failure modes was constructed and the 
corresponding leak rates were determined. The failure modes contributing to a gas leak that were 
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analyzed include either a clean break (flow coefficient (C<*) = 0.98) or a jagged-edged hole (C^ = 
0.68) in the gaseous hydrogen pressurization line and either a clean break (Cj = 0.98) or a jagged- 
edged hole (Cd = 0.68) in the gaseous hydrogen external tank vent. The failure modes contributing 
to a liquid leak that were analyzed include a break in the liquid hydrogen recirculation line and a 
break in the liquid hydrogen feedline, both at or near the external tank interface. Based on this 
analysis, if one were to consider the previously described failures only, a break in the liquid 
hydrogen recirculation line with a larger structural failure occurring in the 72- to 73-second time- 
frame, coinciding with the visibly identified premature ullage pressure decay, seems plausible. The 
analysis was performed in a speculative sense and is not to preclude the possibility that the leak 
could have resulted from a hole or a structural failure elsewhere in the tank. 
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Figure 5 
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.- Average gaseous hydrogen leak rate and relative hole size for each time interval. 
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Relative hole 

diameter, cm (in.) Liquid leak rate, kg/s (Ibm/s) 
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Figure 6.- Average liquid hydrogen leak rate and relative hole size for each time interval. 
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ULLAGE PRESSURE RECONSTRUCTION - DISCUSSION OF ANALYSIS 


The external tank and its pressurization system are described in this section to lay a 
foundation for the detailed discussion of the analysis. With the system description complete, the 
development of the analytical pressurization model used is explained, and then the flight ullage 
pressure reconstruction results are presented. 


PRESSURIZATION SYSTEM DESCRIPTION 

A diagram of the Challenger main propulsion system is shown in figure 7. The system is 
composed of the Orbiter main propulsion system, the Space Shuttle main engines, and the external 
tank. Two separate cryogenic tanks (oxygen and hydrogen) comprise the external tank system. It is 
the hydrogen portion of the external tank and its pressurization system that are of interest. 

A schematic of the hydrogen external tank pressurization system is shown in figure 8. As 
indicated, liquid hydrogen leaves the base of the tank through a 43-centimeter (17 inch) feedline (a) 
and feeds each of three main engines (b). The outlet pressurization gas from each of the three main 
engines passes through a flow control valve (c), where it combines with pressurization gas from the 
other two engines. This combined pressurization gas is channeled into the forward end of the ullage 
and serves to maintain the ullage pressure in a control band of approximately 226 kPa (32.8 psia) to 
230 kPa (33.4 psia). 

The purpose of the flow control valves is to control and limit the flow of pressurization gas 
into the forward end of the ullage. Each flow control valve is represented in the schematic by two 
orifices in parallel. The orifices are always in the choked-flow condition, and each set can simulate 
one of two flow control valve positions. When both orifices are opened, the equivalent flow control 
valve position is full open and permits maximum flow rate. When one orifice is opened and one is 
closed, the equivalent flow control valve position is partly open and limits the flow rate to a 
minimum. 

The flow control valve control system is a closed-loop feedback system with ullage pressure 
being the feedback variable. Each of three ullage pressure transducers (d) acts independently to 
control the three engine outlet flow control valves. When ullage pressure rises above a 
predetermined transducer limit (different for each flow control valve), a signal is sent to close the 
corresponding valve and, consequently, to reduce the pressurization gas flow rate. Likewise, when 
the ullage pressure begins to fall below a predetermined transducer limit, a signal is sent to open 
the corresponding valve and, consequently, to increase the pressurization gas flow rate. The control 
system, therefore, maintains the ullage pressure in a prescribed band, and a plot of ullage pressure 
as a function of time assumes a saw-toothed shape. 
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Orbiter MPS 


Space Shuttle 
main engines 



Figure 7.- Challenger integrated propulsion system. 



Figure 8.- Schematic of hydrogen external tank repressurization system. 
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PRESSURIZATION MODEL DEVELOPMENT 


Overview 


The first objective of this study was to determine whether a leak was, in fact, present. By 
reconstructing hydrogen external tank ullage pressure, a comparison could be made with actual 
measured ullage pressure to rule out a control dynamics malfunction. This reconstruction could be 
used not only to determine whether a leak was present but also to determine its initiation time by 
noting the time at which actual and predicted pressures begin to deviate. The first step in the 
analysis was to create a pressurization model that would enable reconstruction of predicted ullage 
pressure. 


Ullage Pressure Prediction 

The pressure inside the ullage volume is defined by the ideal gas law as 


P(t) = 


T Julk 

V(t) 


( t ) 


( 1 ) 


where 


N(t) = total moles of gas in ullage, kmol (lb-mol) 
Tbulkti) = mean bulk ullage temperature, K (°R) 

V{t) — total ullage volume, m 3 (ft 3 ) 

R = universal gas constant, kN— m 


kmol — K 


/ ft — Ibf \ 
V lb — mol-°R / 


If instantaneous quantities of V, N , R , and are known, predicted ullage pressure can be 

reconstructed. 


Throughout ascent, as liquid hydrogen leaves the tank, pressurization gas expelled from the 
main engines is channeled continuously into the forward end of the ullage in the hydrogen external 
tank. Thus, ullage volume V is a function of the original volume of liquid hydrogen in the tank and 
the amount of liquid hydrogen which leaves the tank. Likewise, the total moles of gas in the ullage 
N is a function of the original number of moles of gas plus the number of moles of pressurization gas 
channeled into the forward ullage. Liquid hydrogen boiloff will add to the total moles of gas; 
however, without knowing the temperature profile within the ullage, the amount of boiloff is 
difficult to predict. Therefore, boiloff was not considered directly but is assumed to be constant 
between flights and inherently compensated for in a bulk temperature correlation. 

Finally, although an ullage temperature measurement exists, this temperature is not 
indicative of a mean bulk temperature. The ullage compartment is somewhat stratified, and the 
temperatures of the relatively warm inlet pressurization gas and the cold liquid hydrogen form the 
bounds on temperature. A measured temperature, therefore, is strongly a function of the transducer 
location and does not indicate the bulk temperature accurately. Consequently, a mean bulk ullage 
temperature T )> u ik, as opposed to a measured temperature, was determined for use in this analysis, 
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and it was derived from a correlation using data from STS 51-L (before 60 seconds) and from 
previous missions. The specifics of the correlation are addressed later in the subsection titled "Mean 
Bulk Ullage Temperature Correlation.” 

The fluid flow diagram of the pressurization system shown in figure 9 indicates the critical 
parameters needed to predict ullage pressure using the ideal gas law of equation (1). Calculated 
parameters were computed as follows: 

1. Gaseous hydrogen mass flow rate m(f), kg/s (lbm/s), assuming choked-flow conditions 


m(i) — 

10.75PC , A 

a 

(SI units) 

Vt 

m(t) = 

0.14 PC J A 

a 

(English units) 

Vt 


(2a) 

(2b) 


where 


SI units 

P = engine outlet pressure, kPa 
T = engine outlet temperature, K 
CdA = effective flow area, m 2 

= 0.3662 cm 2 (engine 1; valve open) 
= 0.0712 cm 2 (engine 1; valve closed) 
= 0.3659 cm 2 (engine 2; valve open) 
= 0.0688 cm 2 (engine 2; valve closed) 


English units 

P = engine outlet pressure, psia 
T = engine outlet temperature, °R 
CrfA = effective flow area, in 2 

= 0.05676 in 2 (engine 1; valve open) 

= 0.01103 in 2 (engine 1; valve closed) 
= 0.05671 in 2 (engine 2; valve open) 
= 0.01066 in 2 (engine 2; valve closed) 
= 0.05636 in 2 (engine 3; valve open) 
= 0.01037 in 2 (engine 3; valve closed) 


= 0.3636 cm 2 (engine 3; valve open) 
= 0.0669 cm 2 (engine 3; valve closed) 


The derivation of equation (2) is given in appendix B. It was assumed that flight data indicated the 
position of each valve accurately. 

2. Total moles of gas N (t) 


N He - ESC + N H 2 -ESC + 


i: 


rh(t)dt 


M*) N He-ESC + N H 2 -ESC + 


i: 


ESC 2.016 kg/kmol 

m(t)dt 

esc 2.016 lbm/lb — mol 


(SI units) 


(3a) 


(English units) (3b) 
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Measured Data 

VP = valve position 
P = engine outlet pressure 
T = engine outlet temperature 
Q □ volumetric flow rate 
P U H = ullage pressure 
T meas = ullage temperature 


Calculated Parameters 

m(t) = gaseous hydrogen 
mass flow rate 
N(t) = total moles of gas 
V(t) = ullage volume 
T bulk = mean bulk temperature 


Figure 9.- Fluid flow schematic of hydrogen external tank repressurization system. 
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where 

N He~ESC = number moles gaseous He at engine start command, kmol (lb-mol) 

Nh 2 -esc = number moles gaseous H 2 at engine start command, kmol (lb-mol) 

m(t) — total gaseous hydrogen mass flow rate as a function of time, kg/s (lbm/s) 

It is assumed that no net condensation and evaporation takes place between the ullage and the 
liquid hydrogen. 

3. Ullage volume V(t) 

™=wC sc ww ' 

where 

V$sc = ullage volume at engine start command, m3 (ft 3 ) 

Q(/) = total engine liquid hydrogen flow rate as a function of time, m 3 /s (ft 3 /s) 

It is assumed that no net condensation and evaporation takes place between the ullage and the 
liquid hydrogen. Appendix C contains initial propellant load data used to determine the initial 
ullage volume prior to main engine ignition. Engine start command occurs 6 seconds before lift-off. 


Initial Reconstruction 


Using the technique just described, an initial reconstruction of predicted ullage pressure as a 
function of time was performed for STS flights 51-L, 61-A, and 51-F. Measured ullage temperature, 
as opposed to a derived mean bulk temperature, was used as the temperature parameter for an 
initial iteration. The purpose of this initial reconstruction was not to produce results but rather to 
provide trends that would later aid correlation of a mean bulk ullage temperature. 

Results of the initial pressure reconstruction are shown in figures 10 to 12. As the figures 
indicate for all three flights, the predicted ullage pressure deviates by 3.4 kPa (0.5 psia) above 
actual ullage pressure at 30 seconds and increases to approximately 13.8 kPa (2.0 psia) above actual 
at 73 seconds. Since pressure is directly proportional to temperature, it follows that a difference in 
actual and predicted ullage pressure parallels the difference in required mean bulk ullage 
temperature and measured ullage temperature. 

To investigate the relationship between the measured temperature and a required mean bulk 
ullage temperature further, a reconstruction of required mean bulk ullage temperature (required 
effective temperature) was performed. Required temperature was computed using actual ullage 
pressure and is equal to the temperature required to match predicted to actual ullage pressures. 
Results are shown in figures 13 to 15. 
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Figure 10.- STS 51-L initial ulla 



Figure 1 1 - STS 61-A initial ulla 
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Figure 14.- STS 61-A comparison of measured and effective bulk temperature. 



Figure 15.- STS 51-F comparison of measured and effective bulk ullage temperature. 
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Two important observations can be made. One, for all three flights, the required temper- 
ature, in a gross sense, follows a measured temperature. This result implies that, perhaps, a mean 
bulk temperature can be correlated with a measured temperature. Second, in all three flights, from 
approximately 30 to 70 seconds, the required effective temperature deviated from approximately 
0 K (0° R) and decreased to 14 K (25° R) below measured temperature. Thus, time seems to affect the 
relationship between measured and bulk temperature as well. With trends of this initial recon- 
struction noted, the analysis proceeded with a focus on correlating a mean bulk temperature such 
that pressure reconstruction for STS flights 61-A and 51-F would coincide with the actual pressure 
history. 

Computer programs used to perform engineering calculations on STS 51-L flight data are 
contained in appendix D. 


Mean Bulk Ullage Temperature Correlation 

The initial reconstruction phase of this analysis revealed that a mean bulk temperature 
required to make predicted pressure equal actual pressure appeared to be a function of both 
measured ullage temperature and time. With this first note, a simple correlation of bulk 
temperature as a function of measured temperature was derived for STS 51-F and reapplied to both 
STS 61-A and STS 51-L. Pressure was then reconstructed using T^ u ik = f ( measured temperature) 
as the temperature parameter in equation (1). This correlation, however, did not prove useful, and 
it was determined that measured temperature alone is an insufficient correlating parameter. 

Since measured temperature alone was not adequate for correlation, the function of time 
needed to be assessed. It is reasonable to postulate that as the ullage volume expands, the location 
of the mean bulk temperature moves farther away from the location of measured temperature such 
that the mean bulk temperature becomes increasingly less than the measured temperature. Since 
ullage volume increases with time, it is not unreasonable to suggest that the function of volume on 
an effective bulk temperature might be manifested, in a rough sense, in the form of time. 
Furthermore, if the throttling profiles, which show the amount of liquid hydrogen leaving the tank 
as a function of time, of all three flights were the same, time alone would be a sufficient correlating 
parameter (for correlations between flights). The throttling profiles are different, however, and at 
any one time, the ullage volume for each of the three flights is not constant. Therefore, volume, not 
time, was chosen as the second correlating parameter. 

With key correlating parameters of measured temperature and ullage volume chosen, a two- 
variable numerical fit for bulk ullage temperature was performed. The goal was to determine a 
single correlation using data from all three flights. The timeframe considered was from 40 to 60 
seconds. The timeframe prior to 40 seconds was not considered since the period of interest was only 
during the final seconds of flight. Data for the period after 60 seconds were not used either since the 
MPS data for flight STS 51-L deviated from nominal during this timeframe. 

The correlation was derived in a two-step fashion. First, a series of least squares fits of 
required bulk ullage temperature (as determined in the initial reconstruction) as a function of 
measured temperature T was determined for constant ullage volume. The measured temperature 
used in these fits was calculated from a separate numerical fit performed to smooth out oscillatory 
flight data for this parameter. The following series of curve fits, one for each ullage volume 
considered, resulted. (They are shown in fig. 16.) The best fit describing the effect of measured 
temperature Ton the bulk temperature T^/* for constant ullage volume was linear. 
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Figure 16 .- Bulk ullage temperature correlation - two-variable Fit (ullage volume constant). 
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V=V u T bulk = f\T) = 

v = v 2 ,T bu[k = f ro = 
v = V 3 , T bulk = f(T) = 

V=V 4 ,T bu i k = f{T) = 

A least squares fit was performed on the coefficients of this series of curves to determine the effect of 
volume. The result was a curve fit for the same order coefficients derived previously, and it is shown 
in figure 17. 

A = f{V) = /fai, 02 , 03 , 04 ) 

B = f(V) = /(61, 62, 63,64) 

Finally, the two-curve fits were combined resulting in a correlation of mean bulk temperature as a 
function of ullage volume and measured temperature. 

T bulk = A + BT, T = K(°R) ( 5 ) 

where 

A = -138.331 + 1.558 X V, V = m3 
B = 2.993 - 0.203 W, V = m3 
A = -248.9958 + 7.955 X 10-2V, V = ft3 
B = 2.993 - 3.4105 X 10-2 W, V = ft3 

This correlation is based on the assumption that ullage volume and measured inlet tempera- 
ture are the significant parameters in correlating a mean bulk temperature. Factors such as accel- 
eration, vibration, and fluid properties, among others, are important heat-transfer parameters as 
well. However, these additional parameters were assumed to be constant between flights and/or to 
be less significant. As demonstrated in the final reconstructions, the accuracy of the derived corre- 
lation, which is based upon measured temperature and ullage volume alone, supports this 
assumption. 
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Final Reconstruction 

With a correlation of mean bulk temperature determined, a final pressure reconstruction was 
performed for each flight. The correlation, derived for flight data from 40 to 60 seconds, was reap- 
plied from 40 to 73 seconds. For the time period of 40 to 60 seconds, the term T in equation (5) was 
calculated from the fit of measured temperature as a function of time used in deriving the correla- 
tion. For the time period of 60 to 73 seconds, a least squares fit of measured temperature was not 
performed and actual measured temperature data were used for T in equation (5). If a curve fit of 
measured temperature had been performed for this latter period, the reconstructed pressure would 
most likely match the actual more closely for STS flights 61-A and 51-F. This procedure, however, 
would have no effect on flight STS 5I-L results since the temperature was constant from 60 to 73 
seconds. 

Results of the reconstructions are shown in figures 18 to 20. Note that for STS flights 51-F 
and 61-A, reconstructed ullage pressure closely matches ±1.7 kPa ( ±0.25 psia), the actual 
throughout the 40- to 73-second period, verifying the analysis. Flight STS 51-L, on the other hand, 
only matches closely from 40 to approximately 64 seconds, at which time a discrepancy between the 
predicted and the actual pressure occurs. 
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Figure 17.- Bulk ullage temperature correlation — two-variable fit (degree of coefficient constant). 



Figure 18.- STS 51-L ullage pressure reconstruction using mean bulk ullage temperature. 
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INVESTIGATION OF THE STS 51-L ULLAGE PRESSURE ANOMALY - 
DISCUSSION OF ANALYSIS 


INTERPRETATION OF PRESSURE DECAY 

An enlargement of the observed anomaly in ullage pressure for flight STS 51-L is shown in 
figure 21. The match between predicted and actual ullage pressure from the time of correlation to 
approximately 64 seconds is consistent with a no-leak assumption. This result implies that the 
actual masses entering and leaving the tank are accounted for and are equal to those calculated 
from measured quantities. However, in the 64- to 73-second timeframe, there is a discrepancy 
between actual and reconstructed pressure. The reconstructed pressure becomes increasingly 
greater than the actual pressure with a maximum deviation of approximately 11.7 kPa (1.7 psia) 
occurring at 73 seconds. 

The fact that the actual ullage pressure was lower than the predicted ullage pressure can be 
explained if either of two phenomena occurred. One argument is that the ullage volume is increased 
faster than is indicated by the volume calculated from measured quantities of liquid hydrogen 
leaving the tank. This explanation is consistent with a liquid hydrogen leak in the tank. Figure 22 
depicts the amount of liquid hydrogen measured to have exited the external tank (based on engine 
flowmeter data) and the amount of liquid predicted to have exited in order to recreate the actual 
pressure. As observed, in the 64- to 73-second timeframe, more liquid would need to exit than 
measured in order to match the actual and predicted ullage pressures. The difference can be 
accounted for by postulating a liquid leak. 

A second possibility is that there was no liquid leak but, rather, that the amount of gas 
pumped back into the ullage was less than measured to have passed through the flow control 
orifices. This theory is consistent with a gas leak either in the pressurization line or in the tank 
itself. Figure 23 indicates the measured amount of hydrogen gas recirculated to the ullage and the 
amount of hydrogen gas predicted to have entered the tank in order to recreate the actual pressure. 
As observed, the net increase in hydrogen gas in the tank predicted, based on the measured 
pressure, is less than the amount of gas that was supposed to have entered based on the flow 
measured to have passed through the flow control orifices along the path to the external tank. The 
difference can be accounted for by postulating a gas leak. 

With the existence of a leak, either gas or liquid, confirmed, the initiation time has to be 
determined. Ullage pressure reconstruction begins to show a deviation between actual and 
predicted in the timeframe of 64 to 66 seconds. The leak initiation time can be established more 
accurately by determining when the pressure discrepancy exceeds the sensitivity of the analysis. 
Here, an assumed gas leak is used to plot the difference between expected and predicted masses (fig. 
24) as a function of time. In the time period considered, this difference oscillates between - 0.23 kg 
( — 0.5 lbm) and + 0.45 kg ( 4- 1.0 lbm). If this range is considered a limit on the sensitivity of the 
analysis, 65.5 seconds is the earliest time at which a leak can be confirmed. 
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MAGNITUDE OF LEAK 


With the presence of a leak and the initiation time determined, the next objective was to 
estimate the magnitude of the leak. This analysis focuses on a comparison of predicted and actual 
pressure rise and decay rates. A plot of measured ullage pressure for one of the three pressure 
transducers during the last several seconds of flight is shown in figure 25. These empirical pressure 
rise and decay rates can, alternatively, be described analytically by taking the partial derivative of 
the ideal gas law with respect to time. 



Differentiating with respect to time. 
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Figure 25.- Example of empirical ullage pressure rise and decay rates 
for transducer number 3. 


Linearizing for each pressure rise and decay, 


where 


A P 




Pav 

= 

average ullage pressure 

v au 

= 

average ullage volume 

Qau 

- 

average LH 2 volumetric flow rate 

AT 


change in bulk ullage temperature 

N »2a, 

- 

total average moles gaseous hydrogen 



average gaseous H 2 mass flow rate 


= 

total moles gaseous helium 

mw„ 2 

= 

molecular weight of hydrogen 


( 13 ) 


Using equation (13), it is possible to predict the ullage pressure rise and decay rates. For a 
no-leak scenario, equality between empirical pressure rise and decay rates and the rates predicted 
from equation (13) would be expected. Physically, this theory infers that the rate of pressurization 
observed should match the rate of pressurization predicted from ullage mass and volume dynamics. 
An inequality arises when the ullage volume and repressurization gas quantities inferred from 
independent flight measurements (such as liquid hydrogen fuel consumption or downstream 
repressurization gas flow rates) and used in equation (13) differ from the actual. This inequality can 
be corrected, however, if one incorporates either a gas leak (m\k denotes mass flow rate) or a liquid 
leak (Qudenotes volumetric flow rate) into the gas and volume terms of equation (13), respectively. 
Solving implicitly for leaks using this method, one can compute the leak rate as shown in the 
following equations. 


Gas Leak (Assumes No Liquid Leak) 


m lk =(A+B + C + D)[N He 
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where 


API 

A = — — actual 
At I 

P Q 



(Note: "i — 1” refers to all time before the timeframe currently considered.) 

Liquid Leak (Assumes No Gas Leak) 

/V +IJ..AI \ 

Q U =(A + B + C + D)(-^ 1 + I V lk J (15) 

where 


AP 

A = - — actual 

At 
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where p is fluid density. 


Results 


Because the transducer sensitivity was limited to ± 0.55 kPa ( ± 0.08 psia), an estimation of 
the empirical pressurization rise and decay rates from the data in figure 25 was made difficult. 
Given, however, that the sampling rate was high enough to indicate the time at which each pressure 
transducer setpoint was exceeded, it was assumed that the midpoints of each vertical step and the 
associated time defined real points of time and pressure. These points were used to calculate 
average pressure rise and decay rates from the start to the end of each of the three time intervals 
considered. In the last timeframe (72.0 to 72.9 seconds), few data points exist to estimate a pressure 
rise and decay rate. In this case, an upper and a lower bound on decay rate were determined and are 
shown in figure 25. Finally, a similar procedure was performed for each of three transducers, and 
the standard deviation o was computed to estimate the error based on the three measurements. 

(The underlying assumption is that there is no error in calculating the rates for each transducer 
using the method just described.) 

Gas and liquid leak rates were computed at an average time for each pressure rise and decay 
slope beginning with the slope from 63.5 to 66.7 seconds. Results are shown in figures 26 and 27 and 
in table 1. If a gas leak is assumed, calculations show an initial leak of 0.28 kg/s (0.62 Ibm/s) ± 41 
percent at 65.5 seconds which increases to between 0.43 and 0.51 kg/s (0.94 and 1.12 Ibm/s) ± 1 
percent in the 72- to 72.9-second time interval. If a liquid leak is assumed, calculations show an 
initial leak of 72.9 kg/s (160.5 lbm/s) ± 41 percent at 65.5 seconds which increases to between 111.6 
and 133.2 kg/s (245.8 and 293.3 lbm/s) ± 1 percent in the 72- to 72.9-second time interval. In both 
cases, the leak remains relatively constant from 65.5 to 72.0 seconds and increases significantly at 
72 to 73 seconds. This sharp increase in leak rates parallels the premature decay that was initially 
identified in the flight ullage pressure data of figure 1. 
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Figure 26.- Average gaseous hydrogen leak rate and relative hole size for each time interval. 
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Figure 27.- Average liquid hydrogen leak rate and relative hole size for each time interval. 
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TABLE 1.- SUMMARY OF LEAK RATES 
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EFFECTIVE FLOW AREA 


If a gas leak is assumed, gas hole sizes can be computed using choked-flow equations. It was 
assumed that a leak occurred in the side of the tank at a location near the mean bulk ullage temper- 
ature and that flow at the hole was single phase. Thus, the effective leak area for a gas leak is 


where 


A = 


m u VT 

Xk bulk 


10.75 PC d 


m \k bulk \ 

0.14 PC d ) 


(17) 


m\k = gas leak rate, kg/s (lbm/s) 

T bulk = mean bulk ullage temperature, K (°R) 

P = ullage pressure, kPa (psia) 

Cd = flow coefficient = 0.68 (jagged hole) 

= 0.98 (clean pipe exit) 

Hole sizes were computed for each leak rate assuming a jagged hole in the side of the tank 
ullage with a flow coefficient of 0.68. These hole sizes are shown in figure 26. An effective hole size 
for this case is approximately 5.6 centimeters (2.2 inches) ± 23 percent in the 65.5- to 72.0-second 
time interval increasing to between 6.8 and 7.4 centimeters (2.7 and 2.9 inches) ± 1 percent in the 
72.0- to 72.9-second time interval. Hole sizes were also computed for a leak through a short stubbed 
tube {Cd = 0.98, clean pipe exit). Finally, to bound these calculations by considering the 
dependence of temperature on the location in the ullage, the term Tbuik in equation (17) was 
replaced with T meas and T sat . Hole sizes using T meas were the same and those using a lower limit of 
T sat resulted in a decrease in diameter of approximately 45 percent. Results of these calculations 
are summarized in table 2. 

If a liquid leak is assumed, liquid hole sizes can be estimated using theory for two-phase 
critical flow through short tubes and orifices (refs. 1 to 3). The exit point was chosen to be in the 
vicinity of the liquid hydrogen feedline at the base of the tank since a structural failure at this 
location is plausible. The azimuthal orientation was determined to be insignificant. A schematic is 
shown in figure 28. 

For a subcooled liquid undergoing two-phase critical flow through a short tube or an orifice, 
the critical flow rate and pressures are described by 


and 
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where 


m t = critical mass flow rate at throat, lbm/s (kg/s) 

v ge = specific volume of vapor at throat equilibrium conditions, ft 3 /lbm (m3/kg) 
v( 0 = specific volume of subcooled liquid at stagnation conditions, ft 3 /lbm (m3/kg) 
s ge - entropy of vapor at throat equilibrium conditions, Btu/lbm-°R (kJ/kg-K) 

S( e = entropy of liquid at throat equilibrium conditions, Btu/lbm-°R (kJ/kg-K) 
dst e = change of entropy of liquid from stagnation to throat equilibrium 
conditions, Btu/lbm-°R (kJ/kg-K) 

dP = change pf pressure from stagnation to throat critical pressure 
A t = throat area, ft 2 (m 2 ) 

N = empirical description of partial phase change occurring at the throat 
x et = quality at throat equilibrium 
P t = throat critical pressure, lbf/ft 2 (N/m 2 ) 

P 0 = entrance stagnation pressure, lbf/ft 2 (N/m 2 ) 


g = 32.2 — (English units only) 

c lbf-s 2 

Equations (18) and (19) can be solved implicitly for P* and A t (or rrt t) given the upstream 
stagnation conditions, P Qt T 0 , and rh t (or A t ). The upstream stagnation pressure was determined 
from the change in pressure head from point 1 to point 2 of figure 28. If one assumes a location in 
the tank as high as point 1, or P 0 = P ullage* the estimated hole sizes decrease by only approximately 
5 percent. The resulting hole sizes are shown graphically in figure 27 and are summarized in table 
2. These hole sizes were computed for a C<* = 0.98. Hole sizes for a = 0.68 are within 6 percent of 
these estimates and are not shown. The exit quality and void fractions are estimated to be 0.03 and 
0.7, respectively, within the accuracy of the model used. (See ref. 2.) However, given the inherent 
uncertainties that still exist in the modeling of two-phase critical flow, the author would like to 
qualify the hole sizes computed as best estimates and specific to the model used. 

Finally, in a speculative sense, a matrix of four possible structural failures was constructed 
and the corresponding leak rates were determined using the same methods and assumptions 
described previously. Results are shown in table 3 and are presented graphically in figures 26 and 
27 for comparisons to leak rates observed. The failures considered are not meant to be inclusive and 
should not preclude the possibility that the leak could have occurred by way of a hole somewhere 
else in the tank. 
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TABLE 2 - SUMMARY OF EFFECTIVE HOLE SIZES 
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Figure 28.- First law of thermodynamics applied to the hydrogen external tank. 
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TABLE 3.- LEAK RATES ASSOCIATED WITH THE POSSIBLE STRUCTURAL FAILURES 


Flight hardware 

Hydrogen leak rates, kg/s (lbm/s) 


Liquid 

Gaseous 


C d = 0.98 

C d = 0.68a 

C d = 0.98b 

LH 2 recirculation 
line — 10-cm (4 in.) diameter 

41.7 (91.8) 

— 

— 

LH 2 feedline — 43-cm (17 in.) 
diameter 

752.8 (1658.2) 

- 

- 

GH 2 pressurization line — 
5-cm (2 in.) diameter^ 

- 

1.3 (2.9) 

1.5 (3.2) 

GH 2 vent line — 13-cm(5in.) 
diameter 

- 

1.5 (3.2) 

2.1 (4.7) 


a Cd = 0.68 corresponds to a jagged hole, 
b Cd = 0.98 corresponds to a clean line break. 
c Assumes no pressurization gas enters for this break. 

As observed, leak rates for the gaseous hydrogen vent line, the gaseous hydrogen 
pressurization line, and the liquid hydrogen feedline well exceed the leak rates observed across all 
time intervals. A break in the liquid hydrogen feedline, however, could account for the leak 
observed in the interval of 65.5 to 72 seconds, especially given the inherent margins for error in 
modeling two-phase flow. A larger structural failure would need to occur in the last time interval 
from 72.0 to 72.9 seconds to account for the large leak rate corresponding to the sharp decrease in 
pressure that was visibly identified as the premature decay in figure 1. Thus, of these four failures 
considered, a scenario of a break in the liquid hydrogen recirculation line at 65.5 seconds with a 
subsequent larger structural failure at 72 to 72.9 seconds seems plausible. 


CONCLUSIONS 


Hydrogen external tank ullage pressure reconstruction for flight STS 51-L shows that a 
deviation between predicted and actual pressure starting at 65.5 seconds is consistent with flight 
data. Furthermore, this pressurization loss can be interpreted as either a liquid hydrogen or a 
gaseous hydrogen leak. If a gas leak is assumed, calculations show an initial leak of 0.28 kg/s (0.62 
lbm/s) ± 41 percent at 65.5 seconds which increases to between 0.43 and 0.51 kg/s(0.94and 1.12 
Ibm/s) ± 1 percent in the 72- to 72.9-second time interval. If a liquid leak is assumed, calculations 
show an initial leak of 72.9 kg/s (160.5 lbm/s) ± 41 percent at 65.5 seconds which increases to 
between 111.6 and 133.2 kg/s (245.8 and 293.3 lbm/s) ± 1 percent in the 72- to 72.9-second time 
interval. In both cases, the leak remains relatively constant from 65.5 to 72.0 seconds and increases 
significantly at 72 to 73 seconds. This sharp increase in leak rates parallels the premature decay 
that was initially identified in flight ullage pressure data. 
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Finally, a comparison of effective hole sizes required to cause the calculated leak rates 
indicates a break in the gaseous hydrogen pressurization line, in the gaseous hydrogen vent line, or 
in the liquid hydrogen feedline is not likely. A break in the liquid hydrogen recirculation line with a 
larger structural failure occurring in the 72- to 72.9-second time interval, the time of the visibly 
identified premature pressure decay, does seem plausible and the most likely of the four modes 
considered. This comparison was performed in a speculative sense and is not meant to preclude the 
possibility that the leak could have occurred somewhere else in the tank. 
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APPENDIX A 

FLIGHT DATA RECONSTRUCTION 


Plots created from the flight data used for flights STS 51-L, STS 61-A, and STS 51-F are 
shown. Tabular data were acquired initially from the NASA George C. Marshall Space Flight 
Center, Slidell Computer Complex. Interpolations to 0.01 second were performed on all data so that 
the files could be meshed at common times to be called upon during the pressure reconstruction. 

One exception is for engine outlet pressure. These plots are not actual flight data and were created 
from logic based on actual engine outlet pressure flight data plots and valve position data. Please 
see appendix D (Program ENGP51L) for the logic used. Table A-l lists the Slidell files accessed. 
Figure A-l shows the corresponding location of each measurement. Figures A-2 to A-43 are plots of 
the data used as a function of mission elapsed time (MET). 


TABLE A-l.- SLIDELL DATA FILES ACCESSED 


No. 

File code 

Description 

1 

E41R1021D 

LH 2 consumption (engine 1) 

2 

E41R2021D 

LH 2 consumption (engine 2) 

3 

E41R3021D 

LH 2 consumption (engine 3) 

4 

V41P1160A 

Engine 1 outlet pressure 

5 

V41P1360A 

Engine 3 outlet pressure 

6 

V41T1161A 

Engine 1 outlet temperature 

7 

V41T1261A 

Engine 2 outlet temperature 

8 

V41T1361A 

Engine 3 outlet temperature 

9 

V41X1661E 

Valve position (engine 1) 

10 

V41X1662E 

Valve position (engine 2) 

11 

V41X1663E 

Valve position (engine 3) 

12 

T41P1700C 

Ullage pressure (1) 

13 

T41P1701C 

Ullage pressure (2) 

14 

T41P1702C 

Ullage pressure (3) 

15 

T41T1705A 

Ullage temperature 
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Figure A-2 - STS 51-L engine 1 fuel consumption (E41R1021D). 
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Figure A-5.- STS 51-L engine 1 outlet pressure (V41P1160A). 
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Figure A-6- STS 51-L engine 2 outlet pressure (created). 
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Figure A-ll.- STS 51-L engine outlet valve 1 position (V41X1661E). 
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Figure A-12 - STS 51-L engine outlet valve 2 position (V41X1662E). 
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Figure A-13 - STS 51-L engine outlet valve 3 position (V41X1663E). 
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Figure A-14.- STS 51-L average ullage pressure (T41P1700C, T41P1701C, T41P1702C). 
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Figure A-22 - STS 61-A engine 1 outlet temperature (V41T1161A). 
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Figure A-25.- STS 61-A engine outlet valve 1 position (V41X1661E). 



-vmnnnnnnnnnnnn n n on / — 1/ ir 




1 1 1 1 1 1 1 r 

0 10 20 30 40 50 60 70 

MET, s 


10.0 

7.5 
5.0 

2.5 
0.0 
- 2.5 
- 5.0 
- 7.5 

- 10.0 

Figure A-26 - STS 61-A engine outlet valve 2 position (V41X1662E). 
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Figure A-27.- STS 61-A engine outlet valve 3 position (V41X1663E). 
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Figure A-28.- STS 61-A average ullage pressure (T41P1700C, T41P1701C, T41P1702C). 
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Figure A-35.- STS 51-F engine 3 outlet pressure (V41P1360A). 
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Figure A-39.- STS 51-F engine outlet valve 1 position (V41X1661E). 
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Figure A-40.- STS 51-F engine outlet valve 2 position (V41X1662E). 
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Figure A-41.- STS 51-F engine outlet valve 3 position (V41X1663E). 


1 




1 



0 10 20 30 40 50 60 70 

MET, s 


234.4 

(34.0) 


231.0 
(33.5) 1 

227.5 
(33.0) -j 



224.1 
(32.5) -\ 


Figure A-42 - STS 51-F average ullage pressure (T41P1700C, T41P1701C, T41P1702C). 
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Figure A-43.- STS 51-F measured ullage temperature (T41T1705A). 
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APPENDIX B 

DERIVATION OF CHOKED-FLOW MASS-FLOW-RATE EQUATIONS 


The total quantity of gaseous hydrogen channeled to the hydrogen external tank ullage was 
determined by integrating the mass flow rate of gas expelled from each of the three main engines. 
No flowmeter data exist for this mass flow rate, however. Instead, the mass flow rate was computed 
using choked-flow equations taken from engine outlet pressure, engine outlet temperature, and 
downstream valve position flight data. The effective flow area, C<*A, was experimentally deter- 
mined for each valve in both the open and the closed positions. The derivation of equation (2), 
gaseous hydrogen mass flow rate, is given herein. 


Main engine outlet 
?tl> Tn 



P*, T*, V *, p* C d A 


Figure B-l.- Choked engine outlet orifice. 


1 . U pstream stagnation pressure and temperature 
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These quantities are determined from engine outlet pressure Pa and engine outlet 
temperature Ta flight data. 

2. Throat temperature T* 


T* 




k + 1 2.404 


= 0.83 


T* = T n X 0.83 (Kor°R) 
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3. Critical throat velocity V* 


4. Throat pressure P* 


5. Throat density p* 


V* = C = VkRT* 




kJ J 
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6. Mass flow rate m 


m = p *V*AC d 

0.1 6P r P = kPa 

= X 69.45 VrTx AC . 

T tl d A 2 

1 n A = m 

01 2 P . P = lbf/ft 2 

= — X169.58VfT X AC ' 

t\. n d . . 2 

A = in 

P t = kPa 
A = m 2 

P t = lbfi'in 2 
A = in 2 


10 75f t l A ^ 
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APPENDIX C 
INITIAL LOAD DATA 


The following data were used as reference material to determine the initial ullage volume 
prior to main engine start. 


TABLE C-l.- HYDROGEN EXTERNAL TANK LOAD DATA 


Parameter 

STS51-L 

STS 61-A 

STS51-F 

External tank total volume 8 
(pressurized), m3 (ft3) 

1508.1 (53153.3) 

1508.1 (53153.3) 

1508.1(53 153.3) 

Liquid hydrogen load at 
engine start command 1 * 

(includes lines and Space Shuttle 
main engines), kg (lbm) 

105 154 (231 617) 

105 342 (232 030) 

105 176 (231 666) 

Orbiter line and Space Shuttle main 
engine loads,* kg (lbm) 

139.4(307.0) 

139.4(307.0) 

139.4 (307.0) 

Propellant bulk density, 3 
kg/m3 (lbm/ft3) 

70.7 (4.42) 

70.7 (4.42) 

70.7 (4.42) 

Ullage volume at engine 
start command, m3 (ft3) 

23.29 (821.0) 

21.30 (750.8) 

22.97 (809.7) 

Moles gas (H 2 and He) at engine start 
command, kmol (lb-mol) 

6.67(14.7) 

6.58(14.5) 

6.22(13.7) 


a Source — Rockwell International Baseline Propellant Inventories (predictions). 
b Source — Marshall Space Flight Center Flight Evaluation Report (reconstructions). 


C-l 













APPENDIX D 

COMPUTER PROGRAMS CREATED 


The following Fortran programs were used to perform engineering calculations on flight data 
for flight STS 51-L. Similar programs were created for flights STS 51-A and STS 51-F as well but 
are not included with this report. Before the programs are presented, they are listed in table D-l. 
The programs used to perform interpolations on tabular flight data and programs used to perform 
minor calculations and correlations also have not been included in this report. 


TABLE D-l.- COMPUTER PROGRAMS CREATED FOR STS 51-L FLIGHT DATA 


Appendix 

Program 

Purpose 

D-l 

ENGP51L 

To create data files of engine outlet 
pressure 

D-2 

M51L 

To calculate gaseous hydrogen mass flow 
rate from each main engine, then 
integrate to determine total mass 
expelled as a function of time 

D-3 

LHV51L 

To integrate liquid hydrogen flow rate 
from the external tank to determine 
total volumetric fuel consumption 

D-4 

51L 

To compute mean bulk ullage temperature 
and predict ullage pressure for 
comparison to measured ullage pressure 

D-5 

LIQLK51L 

To compute hypothetical liquid leaks 
corresponding to pressurization loss 
observed 

D-6 

GASLK51L 

To compute hypothetical gas leaks 
corresponding to pressurization loss 
observed 


Figure D-l illustrates the manner in which flight data were accessed by each of these 
programs. 
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Figure D-l.- Flight data reduction. 
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APPENDIX D-l 
ENGINE OUTLET PRESSURE 


PROGRAM ENGP51L 


This program creates engine outlet pressure profiles for flight 51L. 
Although tabular data for engine outlet pressure exists, this data 
was acquired at 1 Hz only. Since the fluid system itself responds at a 
frequency greater than 1 Hz, interpolations to 0.01 seconds as performed 
on other flight data would be erroneous. Instead, plots of engine 
outlet pressure were used to create logic so that engine outlet pressure 
could be computed as a function of time at 0.01 second intervals. This 
logic is dependent on downstream flow control valve position and was 
created as follows: 


1 ) 


2 ) 


QAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

Caaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

C 
C 

c 
c 
c 
c 
c 
c 
c 

Caaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 

QAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

c 

C DEFINITION OF VARIABLES 

0 *//ViV/V/V.Vjy>VA/ VA/VA'lVA'A' l VyVA'<VA/A/A/ 

c 
c 
c 
c 
c 
c 

CAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 
CAAA A AAA A A AA A A A AAAA A AA AA A A A A A A AA A A AAA A AA AA AAAAAAA AAAAAA AA AA A A A AA AAAA A A A A A A A AA AAA 
C 
C 
C 
C 


First, the plots were divided into like sections and equations 
of lines were written to match plots of flight data assuming no 
change in flow control valve position. 

Second, for each section, pressure was decreased or increased 
in a step-wise fashion for each change in valve position. 

The reasoning is as follows. During a prescribed time 
interval, pressure fluctuates between an upper an lower limit. 
The upper limit is reached when the downstream flow control 
valve is in the closed position, thereby increasing the back 
pressure experienced by the engine. Likewise, the lower limit 
is reached when the downstream flow control valve is in the open 
position, thereby decreasing the back pressure experienced by 
the engine. It was assumed that changes between the upper and 
lower pressure limits were simultaneous with changes in flow 
control valve changes. Since changes in flow control valve 
position can be detected at a frequency >> 1 Hz, this method 
provides a way to more precisely predict changes in outlet 
pressure that would otherwise be missed by the 1 Hz pressure 
measurment itself. 


PA = engine outlet pressure (psia) 

TIME = mission elapsed time (seconds) 

VP = flow control valve position (logical 1 = closed, 
logical 0 = open) 


READ FLOW CONTROL VALVE POSITION DATA 


1 READ <3,*,ERR=998) TIME, VP1 , VP2, VP 3 
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0AAA A A A A A A A A A A A A A A A A A A AAA AAA AAA A AAAAA A A A A* A A AAAAA AAAAA A A AAA AAA A A A At A A A A A A A A A A A A A A 

c 

C MATCH ENGINE 1 OUTLET PRESSURE LOGIC 

C WITH MISSION ELAPSED TIME 

0 «V«W/V<V««/V<y^//W«A<«/V*//W/v//A//v/VVA(/VWA//VA//VA///V/^ 

c 

IF (TIME.LT. -2.0) GO TO 10 

IF ( (TIME. LT. 4.0) .AND. (TIME. GE. -2.0) ) GO TO 20 
IF ( (TIME. LT. 19.0) .AND. (TIME. GE. 4.0) ) GO TO 30 
IF ( (TIME. LT. 22.0) .AND. (TIME. GE. 19.0) ) GO TO 40 
IF ( (TIME. LT. 35. 0). AND. (TIME. GE. 22.0)) GO TO 50 
IF ( (TIME. LT. 39.0) .AND. (TIME. GE. 35.0) ) GO TO 60 
IF ( (TIME. LT. 52.0) .AND. (TIME. GE. 39.0) ) GO TO 70 
IF ( (TIME. LT. 56.0). AND. (TIME. GE. 52.0)) GO TO 80 
IF (TIME. GE. 56.0) GO TO 90 
C 

0AA AAAAA AAAAA A A A A A A A A A A A A A A AAAAA A A A AAAAA AAA A AAA A A A A A A A AAAAA A AAAAA A A A A AAA AAA A A A A A 
CAAAAAAAAAAAAAAA AAAAAAAAAAAAAAAA AAAAAAAAAAAAAAA A A AAA AAAA AAAA AAAAA AAAA AAAAAA AA AAA 


c 

c 

c 


COMPUTE ENGINE 1 OUTLET PRESSURE 




* **** ********* * 


10 PI = 1192 . 86*TIME + 5725.7 
GO TO 999 
20 PI = 3200.0 
GO TO 999 

30 IF (VPl.EQ.l.O) PI = 3340.0 
IF (VP1.EQ.0.0) PI = 3000.0 
GO TO 999 

40 PI = -113.3 * TIME + 5492.7 
GO TO 999 
50 PI = 3000.0 
GO TO 999 

60 PI = -237.5 * TIME + 11312.5 
GO TO 999 
70 PI = 2060.0 
GO TO 999 

80 PI = 322.5 * TIME - 14710.0 
GO TO 999 
90 PI = 3350.0 
GO TO 999 


CAAAA AAAAAAA A A A A AAAAA A A A A A AAAAAAAAAAAAAAAAAAAAAA A A AAAAAAAA AAAAA A A A A A AAA A A A A A A A A A 
C AAA A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A AAA A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A 

c 

C CREATE DATA FILE CONTAINING ENGINE 1 OUTLET PRESSURE VS. MET 

C 

999 WRITE (4,102) TIME,VP1 ,VP2 ,VP3 ,P1 
102 FORMAT ( IX, 4 (F7 . 2 , IX) ,F10 . 2 ) 

GO TO 1 


C~~~ ~~~ rewind tape 4 containing engine 1 outlet pressure and mission 
C elapsed time and prepare for computions of engines 2 and 3 outlet 

C pressure 

C 

998 REWIND 4 
C 

C A AAA A A A A A A A A A A AAA A A AAA A AAAAAAA AAA A AAA AAA A A A A A AAA A AAA AAAAA AAAAAAA A A A A A A A A A A A A A A A 
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0AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

c 

C READ FLOW CONTROL VALVE POSITION DATA 

n >VA//V-VA//VA<'/V'JV<V^' A//VA/^ACA^A///^^A//V/V/V/V A//VACACA/>VA/A^/VA/yV 


2 READ <4,A,ERR=997) TIME , VP1 , VP2 , VP3 ,P1 
C 

CAAAAAAAA AAAAA A A AAA A A A A A AAAAAA A A A A A A A AAAAAA A A AAA A A A A A A A A AAAAA AAAAAAA AAAA A A A A AAAA 
CAA A A A A A A A AAAAA A A AAA A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A AAA A A A A AAAAAA A A A A A A AAA A 

c 
c 
c 


MATCH ENGINES 2 AND 3 OUTLET PRESSURE LOGIC 
WITH MISSION ELAPSED TIMES 






IF (TIME. LT. -2.0) GO TO 100 

IF ( (TIME. LT. 2.0) .AND. (TIME. GE. -2.0) ) GO TO 110 
IF ( (TIME. LT. 19.0) .AND. (TIME. GE. 2.0) ) GO TO 120 
IF ( (TIME. LT. 21. 5) .AND. (TIME. GE. 19.0) ) GO TO 130 
IF ( (TIME. LT. 35.0) .AND. (TIME. GE. 21.5) ) GO TO 140 
IF ( (TIME. LT. 38.0) .AND. (TIME. GE. 35.0) ) GO TO 150 
IF ( (TIME. LT. 51. 5) .AND. (TIME. GE. 38.0) ) GO TO 160 
IF ( (TIME. LT. 57.0) .AND. (TIME. GE. 51.5) ) GO TO 170 
IF (TIME. GE. 57.0) GO TO 180 


£AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 
0A A A A A A A AAA A A AAA A A A AAA A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A AAAAAA A A AAA AAAA AAA A A A A A AAAA A 



COMPUTE ENGINES 2 AND 3 OUTLET PRESSURE 


no engine outlet data exists for engine 2 due to a transducer failure; 
since engine 3 outlet temperature is closest to engine 2 outlet 
temperature, it was assumed that engine 3 outlet pressure logic would be 
closest to engine 2 outlet pressure logic; therefore, engine 3 outlet 
pressure logic was applied to engine 2; computions of engine 2 pressure 
profile call upon engine 2 valve positions, and therefore, make this 
pressure profile unique to engine 2 


100 

P2 

= 1107.14ATIME + 

5314.29 


P3 

= P2 




GO 

TO 996 



110 

P2 

= 3150.0 




P3 

= P2 




GO 

TO 996 



120 

IF 

(VP2.EQ.0.0) 

P2 = 

2900.0 


IF 

(VP3.EQ.0.0) 

P3 = 

2900.0 


IF 

(VP2.EQ.1.0) 

P2 = 

3250.0 


IF 

(VP3.EQ.1.0) 

P3 = 

3250.0 


GO 

TO 996 



130 

P2 

= -120.0* TIME + 

5180.0 


P3 

= P2 




GO 

TO 996 



140 

IF 

(VP2.EQ.1.0) 

P2 = 

2900.0 


IF 

(VP2.EQ.0.0) 

P2 = 

2600.0 


IF 

(VP3.EQ.1.0) 

P3 = 

2900.0 


IF 

(VP3.EQ.0.0) 

P3 = 

2600.0 


GO 

TO 996 



150 

P2 

= -283.3 * TIME + 

12516.7 


P3 

= P2 




GO 

TO 996 
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160 

P2 

= 1750.0 




P3 

= 1750.0 




GO 

TO 996 



170 

P2 

= 209 . 1*TIME 

- ' 

9018.2 


P3 

= P2 




GO 

TO 996 



180 

IF 

(VP2.EQ.1.0) 

P2 

= 3200.0 


IF 

(VP2.EQ.0.0) 

P2 

* 2900.0 


IF 

(VP3.EQ.1.0) 

P3 

= 3200.0 


IF 

(VP3.EQ.0.0) 

P3 

= 2900.0 


GO 

TO 996 




£*A A A AAAA A A AAAAAAAAAAAAAAAAAAAAAAAA AAAAAAAAAA A A AAAAAAAAAAAAAAA A A A A A A A A A A A A A A A A A A 


c 

c 

c 

c 

c 


CREATE DATA FILE CONTAINING ENGINES 1 , 2 AND 3 
OUTLET PRESSURE VS. MISSION ELAPSED TIME 


996 WRITE (7,101) TIME,VP1 ,VP2 ,VP3 ,P1 ,P2 ,P3 
101 FORMAT ( IX , 4 ( F7 . 2 , IX ) , 3 ( FI 0 . 2 , IX ) ) 

C 

QAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 


QAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

c 

GO TO 2 
997 STOP 
END 
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APPENDIX D-2 

REPRESSURIZATION GAS QUANTITY 


This program computes the total mass of hydrogen expelled from 
the three main engines as a function of time. The data file accessed 
contains flight data parameters of: 


1) 

2 ) 

3) 


valve positions (open or closed) for each of three flow 
control orifices 
engine outlet pressures 
engine outlet temperatures 


This data file was created independently from interpolations of flight 
data and includes data points for mission elapsed times of -7 to 73 
seconds at every 1/100 th second. 


PROGRAM M51L 

DIMENSION CDA<3,2>, TMD0T( 10000) 

C 

QA A A A A A A A A A A A A A A A A* AAAAAAAAAAAAAAAAAA A AAAA A A A AAAAAA A A AAAA AAA AAA A AAAA A A A A A A A A A A A A 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

CAA*AAAAAAAAAAAA*AA*AAA**AAAAA*AAAAAA**A**A**AA*A***AA***AAAA*A**At****AA****A**A 
CA A A A A A A A A A AAAA A A A A AAAA A AAAAAAAAAAAAAAAA A AAAAAA AAAAAAAA A A AAAA A A A A A A A A A A A A A A A A A A A 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

qaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 

qaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 

c 
c 
c 
c 


DEFINITION OF VARIABLES 


CDA(x,y) = effective flow control orifice area as given by 

Rockwell International (in2), where x = flow control 
valve (1,2,3) and y = valve position (1 = open, 

2 = closed) 

P# = engine outlet pressure (psia), where # = engine no. 
(1,2,3) 

RMD0T# = gaseous hydrogen mass flow rate (lbm/sec), where 
# - engine no. (1,2,3) 

engine outlet temperature (R) , where # = engine no. 
mission elapsed time (sec) 

integrated mass of gaseous hydrogen expelled from 
the three main engines 

TMDOT = total gaseous hydrogen mass flow rate from the 
three main engines (lbm/sec) 

VP# = flow control valve position (logical 1 = closed, 
logical 0 = open), where # = flow control valve 


SET START CONDITIONS 


T# 

TIME 

TMASS 


CDA (1,1) 
CDA (1,2) 
CDA (2,1) 
CDA (2,2) 
CDA (3,1) 
CDA (3,2) 
TMDOT ( 1 ) 


0.05676 

0.01103 

0.05671 

0.01066 

0.05636 

0.01037 

0.0 
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TMASS =0.0 


0AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

CAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

c 
c 
c 
c 


READ ENGINE OUTLET CHARACTERISTIC DATA 


5 READ (16,*,ERR=999) TIME,VP1 ,VP2 ,VP3 ,P1 ,P2 ,P3 ,T1 ,T2 ,T3 


CAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 
0AAAAA A A A A A A A A A A A A A A A A A A A A A A A A AAAAA A AAAA A A A A A A AAAAAA A A A A A A A A A A A A A A A A A A A A A A A A A A A A 

c 

C ASSIGN EFFECTIVE FLOW AREA FOR EACH SET OF ORIFICES 

C BASED ON CURRENT VALVE POSITION 


IF (VP1.EQ.1.0) J = 2 

IF (VP1.EQ.0.0) J = 1 

IF (VP2.EQ.1.0) K = 2 

IF (VP2.EQ.0.0) K = 1 

IF ( VP3.EQ .1.0) L = 2 

IF (VP3.EQ.0.0) L = 1 

C 

CAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

Caaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 

c 
c 
c 
c 
c 

c 
c 
c 
c 
c 
c 


c 

C" 

c 

c 

c 

c 


COMPUTE GASEOUS HYDROGEN MASS FLOW RATE AND 
INTEGRATE WITH RESPECT TO TIME 


before time = -4.8, it was assumed that no gas was expelled from 
the engines (-4.8 is the average time between a change in engine 
outlet pressure from 0 to > 0 psia; a pressure > 0 indicates that 
gas is being detected at the pressure transducers); the next step 
ensures that no mass flow rate calculations are performed before this 
time 


IF (TIME. LE. -4.8) 
IF ( TIME. LE. -4.8) 
IF (TIME. LE. -4. 8) 


PI = 0.0 
P2 = 0.0 
P3 = 0.0 


gaseous mass flow rate through each flow control valve is computed 
using choked flow equations; the resulting equation for mass flow rate 
is given below and is a function of engine outlet pressure, engine 
outlet temperature and effective flow area only 

RMD0T1 = <0.14*P1*CDA(1,J) >/(Tl**(0.5> ) 

RMD0T2 = ( 0 . 14*P2*CDA( 2 r K) ) / (T2** ( 0 . 5 ) ) 

RMD0T3 = (0.14*P3*CDA(3,L) >/<T3**(0.5> ) 

TMDOT(I) = RMD0T1 + RMD0T2 + RMD0T3 

total mass flow rate is integrated with respect to mission elapsed 
time to determine total mass of gaseous hydrogen expelled 

TMASS = TMASS + 0 . 01*TMDOT( I ) 


QAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

Caaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 
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nnonno 


i 


CREATE DATA FILE CONTAINING TOTAL GASEOUS HYDROGEN MASS 
EXPELLED FROM MAIN ENGINES AS A FUNCTION OF MISSION 
ELAPSED TIME - TO BE ACCESSED BY ULLAGE PRESSURE 
RECONSTRUCTION PROGRAM 




/V^<VA/A/^ 


| 


WRITE (17,101) TIME,RMD0T1 ,RMDOT2 ,RMDOT3 ,TMDOT( I ) r TMASS 
101 FORMAT (1X,6(F10.2,1X) ) 

C 

Q A A A A At A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A Ac A A A A A A A A A A A A A A A At Ac 

c 

GO TO 5 
999 STOP 
END 


i 
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APPENDIX D-3 

LIQUID HYDROGEN CONSUMED 


PROGRAM LHV51L 
DIMENSION TQ( 100000) 

C 

CAAAAAAAAAAAAAAAAA AAAAAAAAAAAAAAAAAAAAAAAA AAAAAAAAAAAAA AAAAAAAAAAAAA A AAAAAAAAA A A 
CAAA A A A A A A A A A A A A A A A A A A A A A A AAAAAAAAAAAAAAAAAAAAAAAAAAA A A A A AAA A A A A A A A A A A A A AAA A A A A A 

c 

C This program computes total volume of liquid hydrogen that leaves 

C the hydrogen external tank and is fed to each of the three main engines. 

C The data files accessed contain the parameters 
C 

C 1) liquid hydrogen fuel consumption (volumetric flow rate) for 

C each main engine 

C 

C These data files were created independently from interpolations on 
C flight data and include data points for mission elapsed times of -7 to 
C 73 seconds at every 1/100 th second. 

C 

CAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 
CAAAAAA A AAAAAAA A A A A A AAAAAAAAAAA A A A A AAAAA A AAAAAA AA A A A AA A A AAAAAAAA A A A A AAA A A A AAAAAA 


C 

C 

c 

c 

c 

c 

C 

C 

C 

C 

c 

c 

c 


DEFINITION OF VARIABLES 


Q# = liquid hydrogen fuel consumption rate per engine (gall/min) , 
where # = 1,2,3 

QFIT = liquid hydrogen fuel consumption rate per engine as computed 
from straight line fit of throttling profile for -5.5 < t 
< -2.0 (gall/min) 

TIME = mission elapsed time (seconds) 

TQ = total liquid hydrogen fuel consumption rate (ft3/min) 

VLH = total liquid volume consumed by all three engines (ft3) 


CAAAAAA A A AAAAAAAA AAAAAAA AAAAAAA AAA AAAAA AAAAAAAA AAAAAA AAAAAAAAAAAAAAAAAAAAAAAAAAA 
CAAAAAA A AAA A AAAAA A A A A A A A A A AAAAA A A AAAAAAAA AAAAAAAAAAA AAA A A AAA A AAA AAAAAA AAAAAAAAA A 


C 

C 

c 

c 


SET START CONDITIONS 


VLH =0.0 
TQ( 1 ) = 0.0 
C 

CAAAAAA A A A A AAA A A A A A A A A A A A A A A AAAAAA AAA A AAA AAAAA AAAAAA AAA A A A A A A A A A A A A A A A A A A A A A A A A A 
CAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 


c 

c 

c 

c 


READ ENGINE FUEL CONSUMPTION DATA 


DO 10 I = 2,100000 
5 READ (29,*,ERR=99) TIME, 01 ,02, Q3 
IF (TIME. GT. 73.0) GO TO 99 


C 

C ~~~~ — throttling profiles (fuel consumption vs. time) show sporadic 
C fluctuations in flow rate prior to t = -2.0 seconds? these 

C fluctuations were attributed to flowmeter start-up transients, 

C therefore, tabular data was not used prior to t = -2 seconds; 
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instead, a straight line fit of fuel consumption as a function of time 
was computed from throttling profiles (-5.5 < t < -2.0) in order to 
wash out transients; this fit is given below 
C 

IF ( (TIME.GE.-5.5) .AND. ( TIME. LE. -2 . 0 ) ) QFIT = 

$ 3 A (4285.71 * TIME + 23971.43) 

IF ( (TIME.GE. -5 . 5 ) .AND. (TIME.LE. -2.0) ) FLAG = 1 
IF (TIME. GT. -2.0) FLAG = 0 
IF ( FLAG . EQ . 1 ) Q1 « QFIT 
IF ( FLAG . EQ . 1 ) Q2 = QFIT 
IF (FLAG.EQ. 1 ) Q3 = QFIT 
TQ ( I ) = <Q1 + Q2 + Q3 ) * 0.002228 
C 

flight data indicates an increase in flow rate above a steady 
C recirculation value (approx. 410 gall/min) for each engine to occur at 

C an average time of t = -5.5 seconds; it was assumed that at this time 

C recirculation of liquid hydrogen ceased and a net flow out of the 

C external tank began; thus, prior to t = -5.5 seconds, fuel consumption 

C is set to zero 

C 

IF (TIME. LT. -5.5) TQ(I) = 0.0 
C 

QAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

QAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

c 

C INTEGRATE FUEL CONSUMPTION WITH RESPECT TO TIME 

C 

integration performed using trapezoidal rule 

C 

IPREV =1-1 

VLH = 0.5*(TQ(I> + TQ(IPREV) )*0.01 + VLH 
C 

CAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 
CAAAAAAAAA A AA AAAAAAAAAAA A A AAAA A AAA A AAAAAA A AA AAA AA AAAAAAAA A A A AA A A AA AAA AAAAA A AAA A A 


c 

C 

c 

c 

c 

c 


CREATE DATA FILE CONTAINING TOTAL VOLUME OF LIQUID 
HYDROGEN CONSUMED FROM EXTERNAL TANK AS A FUNCTION 
OF MISSION ELAPSED TIME - TO BE ACCESSED BY ULLAGE 
PRESSURE RECONSTRUCTION PROGRAM 


C 

WRITE (30,101) TIME, Q1,Q2,Q3, VLH 
101 FORMAT ( IX, F7 . 2 , IX, 4 ( F12 . 3 ) ) 

C 


QAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

CAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 


10 CONTINUE 
99 STOP 
END 


Dll 



APPENDIX D-4 

ULLAGE PRESSURE RECONSTRUCTION 


PROGRAM 51L 


This program reconstructs hydrogen external tank ullage pressure 
for flight STS-51L. The data files accessed contain the critical 
parameters of : 


0AA A AAA A A A A A A A A A A A A A A A A AAA A AAA A AAA AAA A AAA A A A A A A A A A A AAA A A A A A A A A A A A A A A A A A AAA A A A A A A 
0 A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A AAAA A A A AAA AAAAAA A A A A AAA A A A AAAAA 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

C These data files were created independently from interpolations 

C and computations on flight data and include data points for 

C mission elapsed times of -7 to 73 seconds at every 1/100 th sec. 

C 

0A A A A A A A A AAA A A AAA A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A 
0AA A A A A A A A A AAA A AAA A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A AAAA A A A A A A A A AAA A A A A A 


1 ) 

2 ) 

3) 

4) 


volume of liquid hydrogen consumed by main engines 
(computed in program LHV51L) 

mass of gaseous hydrogen measured through flow 
control orifices (computed in program M51L) 
measured ullage temperature 
measured ullage pressure 


C 

REAL L IDEAL, L 

DIMENSION TGLOSS ( 100000 ) , TLL0SS( 100000 ) , SEC(IOOOOO) 

TIME1 '= -7.00 
C 

CAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

CAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 



c 


SET ESC ULLAGE CONDITIONS 




TGLOSS = total gas mass lost (lbm) 

TLLOSS = total liquid volume lost (ft3) 

NESC = total moles of gaseous hydrogen and helium 
at engine start command (moles) 

VESC = ullage volume at engine start command (ft3) 






TGLOSS ( 1 ) =0.0 
TLLOSS ( 1 ) = 0.0 
NESC = 14.7 
VESC = 821.0 


0AA A AAAAAA A A AAAA A AAA A A A AAAAAA A AAA A A AA AAAAAAAAA A A A AAAA A A AAAAAAAAAAAAA A A AAA AAA AAAA 
0AAAAA AAA AAA A A A A A A A A AAA A AAA A AAAAAA A A A A AAA A A A A A A A AAA AAA A A AAAA AAAAAAAAAAA AAAAAA AAA 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


READ ULLAGE CHARACTERISTIC DATA 


TIME# - mission elapsed time (sec), where # =1,2, 3, 4 

V = measured volume of liquid hydrogen consumed from external 
tank (ft3) 

GMASS = mass of hydrogen gas computed to have passed through flow 
control orifices and expected to re-enter the ullage (lbm) 
TMEAS = measured ullage temperature (R) 
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c 

c 

c 

c 


c 


c 


P = measured average ullage pressure (psia) 




‘/A/A/A/A'A'A/A/A/.V/^A 




DO 10 NN = 2,1000 

READ <30,*,ERR=99) TIME1 ,* ,* , * ,V 

READ ( 17 ,* ,ERR=99 ) TIME2,*,*,*,*,GMASS 

READ (43,*,ERR=99) TIME3 ,TMEAS 

READ (39,*,ERR=99) TIME4,*,*,*,P 

IF (TIME1.LT. -5. 5) GO TO 10 
IF (TMEAS.EQ.0.0) TMEAS = 372.0 
IF ( TIME4 . GT .72.9) P = 33.12 
IF (TIME1.GT.73.0) GO TO 99 
IF ( TIME1 . NE . TIME2 ) WRITE (6,100) 

IF ( TIME2 . NE . TIME3 ) WRITE (6,100) 

IF ( TIME3 . NE . TIME4 ) WRITE (6,100) 

100 FORMAT (IX, 'ERROR') 


CAAA A A A AAAAAAAAAAAA AAA A AAAA A A A A A A AAA A A A A A A A A A A A A A A A A AAAA A A A A A A A A A AAA A A A A AAAAA AAA 


CAA AAAAAAAAAAAAAAAAA A A A A AAAA AAA A A A A A A A A A A A A A A A A A A AAA A A A A AAA A A A A A AAAAAAAAAAAA A AAA 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


✓s/ ✓V' A/ /V +sS A/ A/ /V A/ /V A/ /V A/ A/ 


RBAR = 
TVOL = 

TM0LE = 
TBULK = 
PEXP = 


PREDICT ULLAGE PRESSURE 


universal gas constant ( f t-lbf /lbmole-R) 

total ullage volume predicted from original volume at 

ESC and engine fuel consumption data (ft3) 

total moles of gas predicted in ullage based on moles 

at ESC plus quantity of gas measured to have passed 

through flow control orifices (mole) 

mean bulk ullage temperature - computed from correlation 
and is a function of ullage volume and measured 
temperature (R) 

predicted ullage pressure (psia) 


RBAR = 1545.0 

TVOL = VESC + V 

TM0LE = NESC + GMASS/2.016 

IF (TIME1.LT. 40.0) GO TO 56 


C 

C ~~ a curve fit was determined for measured temperature as a function 
C of time (from 40 - 60 sec) in order to smooth measured ullage 

C temperature data for use in derivation of a mean bulk temperature 

C correlation; thus, final computation of a mean bulk temperature using 

C the correlation requires a fitted measured temperature as opposed to 

C an actual measured temperature; this fit is given below 

C 

IF ( (TIME1.GT.40.0) .AND. (TIME1 .LT. 60 . 0 ) ) 

$ TMEAS = 363.401 -.634ATIME1 + . 013289A (TIME1AA2 . 0 ) 

C 

C 

for time greater than 40 seconds, a mean bulk ullage temperature is 
C computed as given by the correlation below 

C 


AA = -248.996 + 7 . 94547E-2*TV0L 
BB = 2.99292 - 3 . 4105E-2A (TVOLAAQ . 5 ) 
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TBULK = AA + BB * TMEAS 
56 PEXP = (TMOLE*TBULK*RBAR)/(TVOL*144.0) 

C 

CAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

G 

c 


COMPUTE TIME HISTORY GAS LOSS 




RNIDEAL= total moles of gaseous hydrogen predicted in ullage 
(computed using measured ullage pressure in ideal 
gas law) (mole) 

TNLOSS = total moles of gas unaccounted for and attributed 
to a leak (mole) 

TGLOSS = total mass of gas unaccounted for and attributed to 
a leak (mass) 

GIDEAL = ideal number of moles of hydrogen that would need to re- 
circulate to ullage to obtain the measured pressure (mole) 






RNIDEAL = (P*144.0*TV0L)/(RBAR*TEMP) 
TNLOSS = TMOLE - RNIDEAL 
TGLOSS (NN) = TNL0SS*2. 016 
GIDEAL = (RNIDEAL - NESC) * 2.016 


C 

QA A A A A A AAA A A A AAAA A A AA A A A A A A A A A A A A A A A A A A AAAAAAAAA A A A A A AAAAAAAAAAAAAAAAAAAAA AAAA A A 
CAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 




COMPUTE TIME HISTORY LIQUID LOSS 


VIDEAL = predicted ullage volume (computed using measured ullage 
pressure in the ideal gas law) (psia) 

TVLOSS = total volume of liquid hydrogen unaccounted for and 
attributed to a leak (ft3) 

TLLOSS = total mass of liquid hydrogen unaccounted for and 
attributed to a leak (ft3) 

LIDEAL = ideal volume of liquid hydrogen that would need to 
have left the external tank to obtain the measured 
pressure (ft3) 

L = ideal mass of liquid hydrogen that would need to 
have left the external tank to obtain the measured 
pressure (lbm) 


VIDEAL = ( TBULK *RBAR*TMOLE ) / ( P* 144 . 0 ) 
TVLOSS = VIDEAL - TVOL 
TLLOSS (NN) = TVLOSS *4. 42 
LIDEAL = (VIDEAL - VESC) * 4.42 
L = V * 4.42 


C 

0AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 


c 

c COMPUTE EFFECTIVE TEMP FOR NO LEAK 

C 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


A//«/A/A//V/V/V/V/V/sfAC/V/WA/A'A 


w'a^/v'a/a/a^wa/a/a/a//v<v/va//v<v>vA 


TEFF = effective bulk ullage temperature required to make 
measured and predicted ullage pressures match (this 
parameter was used in deriving a correlation for 
mean bulk ullage temperature by setting TEFF = TBULK = 
f ( TMEAS , TVOL ) ) 








TEFF = (P*144.0*TVOL)/(TMOLE*RBAR) 


C 

Q A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A AAA ^ A A A A A * A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A 
0AA A A A AAAAA AAAAA A A A AAAA A A A A A A A A AAAAAAAAA AAAAAAAAAAAAAAAA AAA A AAAAAAAA AAAAAAA A A A A A 


c 

C CREATE PLOT FILES 

C 

sends data every 0.05 seconds to plot files 
C 

A = NN 
B = A/5 

ICHECK = ANINT(B*100.0) 

I = (NN/5) * 100.0 
IF ( ICHECK. NE. I) GO TO 10 
WRITE (84,105) SEC(NN), TEMP 
WRITE (86,105) SEC(NN) , TEFF 
WRITE (87,105) SEC(NN), P 
WRITE (88,105) SEC(NN), PEXP 
WRITE (93,105) SEC(NN) , LIDEAL 
WRITE (90,105) SEC(NN) , L 
WRITE (91,105) SEC(NN), GIDEAL 
WRITE (92,105) SEC(NN), GMASS 
WRITE (82,105) SEC(NN), TGLOSS(NN) 

WRITE (83,105) SEC(NN) , TLLOSS(NN) 

105 FORMAT (IX, F7 . 2 , 1X,F10 . 2 ) 

C 

0 A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A: A A A A A A A A A A A A A A A A A A A A A A A A A 
QAA A A A A A A A A A A A A A A A A A A A A A* A A A A A AAA A A A A A A AAA A AAAAA A A A A A AAAAAAA A A A A* A A A A* A A A AAA AAA A 

c 

10 CONTINUE 
99 STOP 
END 
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APPENDIX D-5 

LIQUID HYDROGEN LEAK RATE 


PROGRAM LIQLK51L 
C 

C*AAAA A AAAA A A AAA A AAAAAAA AAAAAAA A AAA AAA A A A AAA A A A AAA A AAA AAAAAA A AAA AAA A A A A A A A AAA A A A 
CAA A AAAAA A A A A A A A A A A A A A A A A AAAAAAAAAAAAAAAAAAAAAAAAAAAAA A A A A A A A A A AAA A A A A A A A A A A A AAA 

c 

C This program computes hypothetical liquid leaks (assuming 

C no gas leak) that are required to make empirical ullage pressure rise 

C and decay rates match predicted values. Time intervals correspond to 

C either a single decay or rise slope or part of a single decay or rise 

C slope on the saw-toothed ullage pressure plots. The data file accessed 

C contains the critical parameters of : 

G 

C 1) time interval considered 

C 2) empirical ullage pressure rise and decay rates 

C 3) predicted ullage pressure assuming no leak 

C 4) mean bulk ullage temperature 

C 5) rate of change of mean bulk ullage temperature 

C 6) average total mass of hydrogen fed to the external tank 

C ullage 

C 7) mass flow rate of gaseous hydrogen from engines to ullage 

C 8) volumetric flow rate of liquid hydrogen to main engines 

C 9) average total volume of liquid hydrogen consumed by 

C main engines 

C 


CA A AAA A A A A A A AAA A A A A AAAAA A A A A A AAAAAAAAAAAAAAAAAAAAAAAAAAA A A AAA AAA A A AAA A A A A A A AAAAA 
C AAA A A A AAA A A A A AAA A AAAAAAA A A A AAA AAAAAAAAAAAAAAAAAAAAAAAAAAAAA A A A A A A AAAAA A A A AAA A A A 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


DEFINITION OF VARIABLES 


ANSW = computed average liquid leak rate during time interval 
considered (ft3/sec) 

DELT = magnitude of time interval (seconds) 

P = average ullage pressure during time interval 
considered (psia) 

PSLOPE = empirical ullage pressure rise and decay rate (as 

determined from the slope of a straight line fit of 
ullage pressure over the time interval considered) 

( psia/sec ) 

Q = volumetric flow rate of liquid hydrogen to the 
main engines (ft3/sec) 

QLK = assumed liquid hydrogen leak rate (lbm/sec) 

RMASS = average total mass of hydrogen fed to the external 
tank ullage 

RMDOT = gaseous hydrogen mass flow rate to ullage (lbm/sec) 

T = average mean bulk ullage temperature during time 
interval considered 

TSLOPE = rate of change of mean bulk ullage temperature 

(as determined from the slope of a straight line 

fit of computed ullage temperature over the time interval 

considered) (R/sec) 

V = average total volume of liquid hydrogen consumed by 
the main engines (ft3) 

VPOT = total volume of estimated liquid leakage prior to 
time interval considered (ft3) 
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c 

CA AAA AAA A AAAAAAAA A AAAAAAAAAA A A A A A A A A A A A A A A A A A A A A A A A A A A A A A A AAAAAAAAA AAAAAAA A A A A A A 
C AA AA AAA AAA AAA AAA A A A AAA A A A A A A A AAA A A AAA A A A A A A A A A A A A A A A A A A A A AAA A A A A A A A A A A A A A A A A A A A 

c 

C READ INPUT DATA 

C 

READ (49,*,ERR=50)DELT,PSLOPE,P,T,TSLOPE,RMASS,RMDOT,Q,V,QLK,VPOT 
C 

CAAAAAAAAAAAAAA AAAA A A AAAAAAA A A A A A A A A A A A A A A A AAA A A A A A A A A A A A A AAA AAAAAAAAAAAA A A AAA A A 
QAA AAAAA A A AAAAAA A A A AAAA A A A A A A A A A A A A A A A A A A A A AAA A A A A A A A A A A A A A A A AAAAAAAAAAA A A A A A A A A 

c 

C CALCULATE LEAK 

c 

C~ ~ the equation used to compute the liquid leak is not explicit for the 
C leak itself; therefore, iterations were performed by assuming a leak 

C ( QLK ) , computing a leak (ANSW) which is dependent on QLK and then 

C comparing the leak assumed (QLK) to the leak computed (ANSW); the 

C procedure was repeated until QLK = ANSW 

C 

A = PSLOPE A 144.0 
B = (P*TSLOPE*144.0)/T 
E = 14.7 + RMASS/2 . 016 
C = <P*RMDOT*144.0)/(EA2.016) 

U = V + QLKADELT/2.0 + VPOT + 821.0 
D = (PAQA144.0) /U 

ANSW = (( -A + B + C - D)*U)/(P*144.0) 

WRITE (6,99) A,B,C ,D,E,ANSW 
99 FORMAT ( IX, 5 ( F9 . 3 , IX) ,E12 . 4 ) 

50 STOP 
END 
C 

CAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 
CAAAAA A A A AAAAA AAAAAAAAAAA A AAAAA AAA A AAA A A AAAA A AAAA AAAA AAAAAAA AAAAAAAAAAAA A A AAAAAA 
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APPENDIX D-6 

GASEOUS HYDROGEN LEAK RATE 


PROGRAM GASLK51L 


This program computes hypothetical gas leaks (assuming no liq- 
uid leak) that are required to make empirical ullage pressure rise and 
decay rates match predicted values. Time intervals correspond to either 
a single decay or rise slope or part of a single decay or rise slope on 
the saw-toothed ullage pressure plots. The data file accessed contains 
critical parameters of : 


the 


CAAAA****AA**AAA*AAAAAAA*AAAAAAAAAA*AAAAAA*A**A* ************ A* ****************** 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


1 ) 

2 ) 

3) 

4) 

5) 

6 ) 

7) 

8 ) 
9 ) 


time interval considered 

empirical ullage pressure rise and decay rates 

predicted ullage pressure assuming no leak 

mean bulk ullage temperature 

rate of change of mean bulk ullage temperature 

average total mass of hydrogen fed to the external tank 

ullage 

mass flow rate of gaseous hydrogen from engines to ullage 
volumetric flow rate of liquid hydrogen to main engines 
average total volume of liquid hydrogen consumed by 
main engines 


C***AA*AAAAAA*AA*AAAAAAAAAAAA**AA***AA*AA*AAAAA*AAAAA*AA*AAAAA*AAAAA*AAA*AAAAAA 

C***aaaaaaaaaaaaaaaaaaaaaaaa****aaa***************************aaaaaaaaaaaaaaaaa 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


DEFINITION OF VARIABLES 


ANSW * computed average gas leak rate during time interval 
considered (ft3/sec) 

DELT = magnitude of time Interval considered (seconds) 

P = average ullage pressure during time Interval 
considered (psla) 

PLSOPE = empirical ullage pressure rise and decay rate 

(as determined from the slope of a straight line 
fit of ullage pressure over the time interval 
considered ) ( psia/ sec ) 

Q = volumetric flow rate of liquid hydrogen to the 
main engines (ft3/sec) 

RLK = assumed gaseous hydrogen leak rate (lbm/sec) 

RMASS = average total mass of hydrogen fed to the external 
tank ullage (lbm) 

RMDOT = gaseous hydrogen mass flow rate to ullage (lbm/sec) 
RPOT = total volume of estimated gas leakage prior to 
time interval considered (ft3) 

T = average mean bulk ullage temperature during time 
interval considered (R) 

TSLOPE * rate of change of mean bulk ullage temperature 

(as determined from the slope of a straight line fit 
computed bulk ullage temperature over time interval 
considered) (R/sec) 

V = average total volume of liquid hydrogen consumed by 
the main engines (ft3) 
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c 

QAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

CAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 

c 

C READ INPUT DATA 

£ VW/V/VA»^/V^/V<VA»A//VA/A/ 

c 

READ <49,*,ERR=50)DELT,PSLOPE,P,T,TSLOPE,RMASS,RMDOT,Q,V,RLK,RPOT 
C 

C^aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 

CAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 


c 

c 

c 

c 

QA'A'/VA'/VA/ 

c 

c 

c 

c 

c 


CALCULATE LEAK 


the equation used to compute the gas leak is not explicit for the leak 
itself; therefore, iterations were performed by assuming a leak ( RLK ) , 
computing a leak (ANSW) which is dependent on RLK and then comparing 
the leak assumed (RLK) to the leak computed (ANSW); the procedure was 
repeated until RLK = ANSW 


A = PSLOPE a 144.0 
B = (P A Q A 144 . 0 ) / ( V + 821.0) 

C = (PATSL0PEA144.0)/T 

E = ( (RMASS - RPOT - (RLKADELT) /2 . 0 ) /2. 016 + 14.7)A2.G16 
D = (PARMDOTA144.0)/E 

ANSW = ( ( -A - B + C + D) * E) / ( P A 144.0) 

WRITE (6,99) A,B,C,D,E,ANSW 
99 FORMAT ( 1X,5(F9 . 3 , IX) ,E12 . 4 ) 

50 STOP 
END 


CAAAAA A AAAAAAAAAAAAAAAAAAAAAAAAAAAAA AAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 
0A A A A A A AAAA A A A A A AAA A A A A A A A A A A A AAA A A A A A A A A A A A A A A AAAA A A A AAAAA A A A A A A A A A AAAAA A A A A A A 
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